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Abstract. The Boltzmann collision operator for a dilute granular gas of inelastic rough hard spheres is much more intricate 
than its counterpart for inelastic smooth spheres. Now the one-body distribution function depends not only on the translational 
velocity v of the center of mass but also on the angular velocity (0 of the particle. Moreover, the collision rules couple v 
and (0, involving not only the coefficient of normal restitution a but also the coefficient of tangential restitution j3. The 
«/j , aim of this paper is to propose an extension to inelastic rough particles of a Bhatnagar-Gross-Krook-like kinetic model 

Cn ■ previously proposed for inelastic smooth particles. The Boltzmann collision operator is replaced by the sum of three terms 

representing: (i) the relaxation to a two-temperature local equilibrium distribution, (ii) the action of a nonconservative drag 

^ force F proportional to v — u (u being the flow velocity), and (iii) the action of a nonconservative torque M equal to a linear 

^4— ( , combination of (0 and ii (Q, being the mean angular velocity). The three coefficients in F and M are fixed to reproduce the 

O ■ Boltzmann collisional rates of change of Ci and of the two granular temperatures (translational and rotational). A simpler 

version of the model is also constructed in the form of two coupled kinetic equations for the translational and rotational 
velocity distributions. The kinetic model is applied to the simple shear flow steady state and the combined influence of a and 
j3 on the shear and normal stresses and on the translational velocity distribution function is analyzed. 
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o 

^ . INTRODUCTION AND MOTIVATION 

fSJ ■ As is well known, the master equation for the dynamics of an ordinary gas in the rarefied regime is the Boltzmann 

^ ' equation [1]. In the particular case of hard spheres, it reads 

o ■ 

p. drf{r,y;t)+yVf{r,y;t)^Ay\f{r,-,t)J{r,-,t)], (1) 

!>■ ; 

^^ I where /(r,v;f) is the one -particle velocity distribution function and 

o 
o 



r'[yi\f{r,-u)J{r,-u)]^o^Jdy2jdd&{g-d){g-d)[f{ry(-t)f{ryi-t)-f{r,yi-,t)f{r,y2-,t)] (2) 



is the Boltzmann collision operator. Here, a is the diameter of a sphere, &{x) is Heaviside's step function, CT is a unit 
vector directed along the centers of the two colliding particles, g = Vi — V2 is the relative velocity, and the double 
S^ i primes on the velocities denote the initial values {v'/,V2} that lead to {vi,V2} following an elastic binary collision: 

.P.'. v'/ = vi-(a-g)a, v'2'-v2 + ((y-g)<T. (3) 

The basic properties of the Boltzmann colhsion operator are 

'fl'v7'='[v|/,/]=0, fdyyj'\y\fj]=0, (4) 



dyy^r'[y\fJ]^Q, (5) 

3/2 / my2 



■/''[v|/o,/o]=0^/o(v)=«(^^) expf-^^j, (6) 



where V = v — u is the peculiar velocity. Henceforth, for the sake of notation simplicity, the arguments r and t will 
be implicitly understood and generally omitted. Equations (4) and (5) express the collisional conservation of mass. 



momentum, and energy, while Eq. (6) indicates that collisions do not change the local equilibrium distribution function. 
In the latter quantity, m is the mass of a particle and n = Jd\f{\), u = n^' Jd\\f{\), and T = {m/3n) J d\V^f{\) 
are the local number density, flow velocity, and temperature, respectively. Note that, as usually done in the literature on 
granular gases, the Boltzmann constant has been absorbed in the definition of temperature, so that the latter quantity 
has dimensions of energy. 

The mathematical intricacy of the Boltzmann operator (2) makes it difficult to get explicit nonequilibrium solutions 
of Eq. ( 1 ), especially beyond the scope of the Chapman-Enskog method. This has motivated the proposal of simplified 
model kinetic equations [1] where J''^[y\f,f] is replaced by a more tractable term 7r'''[v|/] that preserves Eqs. (4)-(6). 
The simplest and most celebrated of these model kinetic equations is perhaps the one proposed in 1954 by Bhatnagar, 
Gross, and Krook (BGK) and, independently, by Welander [2]. It consists of the replacement 

7«i[v|/,/]^/:^'[v|/]=-v[/(v)-/o(v)], V = y CT2„y^^7>, (7) 

where v(r,f) is an effective velocity-independent collision frequency. Its expression is not fixed by Eqs. (4)-(6) 
and thus it can be freely chosen to optimize further agreement with the original Boltzmann equation. The explicit 
expression in Eq. (7) guarantees that the Navier-Stokes shear viscosity obtained from the BGK model coincides with 
that obtained from the Boltzmann equation in the first Sonine approximation [1]. The BGK kinetic model has proven 
to be much more reliable than its formulation might anticipate, even in far from equilibrium states [3]. 

Needless to say, the complexity of the Boltzmann equation increases when it describes a granular gas made of 
inelastic hard spheres [4]. In that case, the Boltzmann collision operator is given by 

-/'"'^'[vil/,/] = a^jdy2Jdd&{g-a) (g-a) [a-V(v'i')/(v'2')-/(vi)/(v2)] , (8) 

where a < 1 is the coefficient of normal restitution (here assumed to be constant) and Eq. (3) must be replaced by 

v'/ = vi :^ — {o-g)a, V2=V2 + — ^^ — {o-g)a. (9) 

Mass and momentum are still conserved, so that Eq. (4) is also satisfied by 7'"'''. On the other hand, the kinetic energy 
is dissipated by collisions, so that one has 

C^-^/^vv2/-'[v|/,/]>0 (10) 

instead of Eq. (5). Equation (10) defines the so-called cooling rate ^ as a functional of /. Moreover, inelasticity 
prevents the system from reaching an equilibrium state and so Eq. (6) is not fulfilled by 7'"'^'. Instead, an undriven 
homogeneous system freely cools down according io d^T ~ —C,hT, where ^/, is the cooling rate corresponding to the 
so-called homogeneous cooling state //j, which is a similarity solution of the Boltzmann equation depending on time 
only through the temperature. As a consequence, the role of /o is played by //, and Eq. (6) is replaced by [5] 

-/'"^'M//„//,] = ^|--(v//,). (11) 

2 ay 

As in the elastic case, it is desirable to construct a kinetic model for granular gases by replacing the true operator 
y'"'''[v|/,/] with a simpler term /r'"'^'[v|/]. Inspired by the BGK model (7) and by the properties (10) and (11), the 
following kinetic model was proposed in Ref. [5]: 

■/'"^'[v|/,/] ->^'"'='[v|/] - -A(a)v[/(y) -/„(y)] + ||- . [V/(y)] , (12) 

2 ay 

where //, is here the local form of the homogeneous cooling state solution, obtained from the latter by replacing v ^ V 
and by replacing the temperature and density with their local values for the nonequilibrium state under consideration. 
As also happened in the original BGK model (7), the modified collision frequency Xv remains a free parameter. By 
keeping v as given by Eq. (7), the factor X{a) can be chosen, for instance, to optimize agreement with the viscosity, 
the thermal conductivity, or the self-diffusion coefficient for the Boltzmann equation [5, 6]. The latter criterion simply 
yields X{a) = (1 + a)/2. The kinetic model (12) is still too complicated because it requires, on the one hand, the 



solution of Eq. (11) and, on the other hand, the cooling rate as a functional of / through Eq. (10). In the spirit of a 
simple kinetic model both conditions can be relaxed [5] by substituting in Eq. (12) the unknown local distribution 
function //j with the local equilibrium distribution /o and the detailed cooling rate C, with its local equilibrium estimate 
C^^(l-a2)v. 

While describing a granular gas, both the Boltzmann operator (8) and its associated BGK-like model (12) assume 
that the inelastic hard spheres are smooth. However, there exist a number of relevant effects in granular gases where 
roughness plays a crucial role [7, 8, 9]. The main purpose of this paper is to extend the BGK-like model (12) to the case 
of inelastic rough hard spheres characterized, apart from a, by a coefficient of tangential restitution j3. The situation 
is now much more complex since, in addition to the translational degrees of freedom, there also exist rotational ones 
[8, 9]. In fact, one can define two independent granular temperatures, T^ and T™', and thus there are two associated 
energy production rates, t,"^ and <^™', instead of a single cooling rate C,. The proposed kinetic model will be further 
simplified by focusing on the translational distribution function. This latter model will be applied and solved to the 
steady simple or uniform shear flow. 

THE BOLTZMANN EQUATION FOR INELASTIC ROUGH HARD SPHERES 

Let us consider a granular gas made of inelastic rough hard spheres characterized by a constant coefficient of normal 
restitution a and a constant coefficient of tangential restitution j3 . The coefficient a ranges from a = (collisions 
perfectly inelastic) to a = 1 (collisions perfectly elastic), while the coefficient j3 ranges from j3 == — 1 (spheres perfectly 
smooth) to j3 = 1 (spheres perfectly rough). Except if a = 1 and j3 = ±1, energy is dissipated upon collisions. The 
one-particle distribution function /(r, v, O; f ) depends not only on the translational velocity v of the center of mass but 
also on the angular velocity fi). In the dilute regime / obeys the Boltzmann equation 

5,/(r, V, a);f) + V • V/(r, V, fi);f) = /""gi^iv, fi)|/(r, •, sOJlr, ', sO], (13) 

where the collision operator is 

r"^''[x,,(0,\fj]^a^ j dy2J d(02J dd@{g-a){g-d)[{afi)-^f{yl(o'i)f{yl^^^^ 

(14) 
The restituting collision rules are [8, 10] 



where 



2 2 

v7 = Vi-C, V2=V2 + C, (Bi=(Bi CTXC, (Bo = ©2 & X C, 

OK OK 


4/ ^ \ + a-\^ .^ K l+j3-i 
.a- ^= 2 (^•^)^+l + . 2 


g-{a-g)a-^dx{a)i + a>2) 



(15) 



(16) 



/ being the moment of inertia. The value of K depends on the mass distribution within the sphere and runs from the 
extreme values k — (mass concentrated on the center) to K" = | (mass concentrated on the surface); if the mass is 
uniformly distributed, then K" = |. Mass and linear momentum are conserved by collisions, so one has 

|rfv|rffi)/™s'^[v,fi)|/,/]=0, ldyldmyr-^''[y,m\fJ]=0. (17) 

On the other hand, except in the perfectly smooth case (j3 = — 1), collisions tend to decrease the mean angular velocity 
i2(r,f ). This can be characterized by a "de-spinning" rate i^n defined by 

dy [do)a)J"'"^^[y,a)\f,f] = -^anCl, Q.= -[dy f da)a)fiy,a)). 



(18) 

In general, neither the translational nor the rotational kinetic energies are conserved by collisions. This can be 
characterized by the partial energy production rates 

^^ ^-^^ j dy j dOiy^r^^^yMfJl ^'"' ^-.^^ j dy j daico^^^^^yMfJ]- (19) 



Here, T"^ = {m/3n) J d\ J daiV^f{\,(0) and T™' = {I /3n) J d\ J d(0(0^f{\,(0) are the translational and rotational 
temperatures, respectively. In general, collisions produce a transfer of energy between the translational and rotational 
degrees of freedom and, consequently, the energy production rates £, "^ and £, ™' do not have a definite sign. On the other 
hand, unless a = 1 and j3 = ±1, the total energy (translational plus rotational) is dissipated and thus the net cooling 
rate 

C = ^(7'"'<?"' + 7^"'<?"% T = ^, (20) 

is positive definite, where T is the total temperature. 

Note that, instead of T™', one could have alternatively adopted T™ = {I/3n) Jd\ Jd(0 {(0 — Q.) f{\,(0) — 
r™' (1 — X), with X = Kma^Q? / 12T'^°\ as the definition of the rotational temperature. Its associated production rate 

is <^ = (<^''°' — 2^nX)/(l —X). However, the disadvantage of this alternative choice is that, in contrast to ^, the 
"cooling" rate ^ = {^T — 2^qX) /T associated with the alternative total temperature T = ^(r"' + T™ ) = T — ^T^°^X 
is not positive definite and in fact becomes negative in the perfectly elastic and rough case (a = 1, j3 = 1). 

Equations (18)-(20) define the de-spinning rate (Ca)^ the energy production rates (^"^, t,™^), and the cooling rate (^) 
as functionals of the velocity distribution function /. They can be estimated in terms of the local values of n, T"^, T™', 
and Si by the replacements 

/(y,a)) ^r(y)r'(fi)), f{y)^f^{y), 

where 

f(y) = jd(Of(y,(0), r\(0) =^-jdyf(y,(0) 

are marginal distribution functions and f^ is given by Eq. (6) with T -^ T"'. The results are [8, 10] 



(21) 

(22) 
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r°'(i+x; 



'T'tr 
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'T'tr 



(23) 
(24) 
(25) 
(26) 



1 + K- v^ r"' 

Before closing this section it is worthwhile noting that the vanishing of the second integral in (17) is not verified 
in the case of the operator /f™^ [/i 1/2] describing the collisions of tagged particles (label 1) with untagged particles 
(label 2). Under conditions milder than Eq. (21) one has [10] 



1 + K- 2 

(27) 

According to Eq. (27), the average collisional transfer of momentum from component 2 to component 1 in the inelastic 

rough case is A(a, j3) times the value in the elastic smooth case. This property is relevant for self-diffusion problems. 



KINETIC MODELING 



Joint distribution function 



As in the smooth-sphere case, both elastic [cf. Eq. (7)] and inelastic [cf. Eq. (12)], the idea behind a kinetic model 
is the replacement of the complex Boltzmann collision operator y^s^ \j ^ y j \y^ ^ much simpler term 7r™"8'^ [/, /] that 
otherwise retains a number of basic physical conditions. As a natural extension of Eq. (12), the BGK-like model 
proposed here for the joint distribution / is 



r^^\yMfJ\^^'''''^\yMf\ = 



<^n- d 



-A(a,/3)v[/(y,a))-/5'(v)/o™'(a))]+^^.[V/(y,a)) 

/(v,fi))}, 



■|[2Caft + <? (fi)- 



2 da 



Q.) 



(28) 



where 



foiy)^n 



2nT^' 



3/2 



exp 



2V 



^"'""-liS-)'""' 



2r°' 



(29) 



are the local equilibrium distributions at independent temperatures. The second term on the right-hand side of Eq. (28) 
can be interpreted as representing the action of an external drag force F = —j^^Y. Likewise, the third term represents 



an external torque of the form M 



. By construction, the model (28) complies with the 



i[2Cnii + <^ (eo-ii) 

exact properties (17)-(19), regardless of the choice of the effective collision frequency Xv. Based on the property 
(27), here the adopted choice for v and X is given by Eqs. (23) and (27), respectively. As for ^a, S,^, and <^™', they are 
explicitly given by Eqs. (23)-(25). 



Marginal distributions 

The kinetic model (28) represents a significant simplification with respect to the original Boltzmann collision 
operator (14). However, it is still too complicated a model if one wants to apply it to situations where the BGK- 
like model (12) for inelastic smooth spheres allows for exact solutions [6, 11, 12]. In order to get explicit results that 
might be useful to assess the influence of roughness on the basic physical properties, it seems desirable to use Eq. (28) 
as the starting point for an even simpler kinetic model. 

More explicitly, we consider now the kinetic equations for the marginal distributions defined by Eq. (22). When the 
replacement (28) is inserted into Eq. (13) and integration over 0) or over v are carried out, one obtains 



^r/^W+V Vr(v) = -Xia,j3)v [f(y)-f^iy)]+^±.[Yfiy)] , 

1 d 



(30) 



d,r\(o)+u-vr\(o) = -A(a,j3)v [r'{(o)-fr{(o)] 



^^^ {[2CaCl + C\(0-Cl)]r\(0)}. (31) 
Upon writing Eq. (31) we have taken into account that dtti + V • (nu) = and have introduced the approximation 



d\\fi\,0))^nuf°'{0)). 



(32) 



Thanks to this approximation, the kinetic model (28) for the joint distribution / yields the set of two coupled kinetic 
equations (30) and (31). 

Actually, the kinetic equation (30) for the translational distribution function /'' is coupled, via S,^, to the rotational 
distribution /™' only through the first and second moments CI and T™'. Therefore, only the equations for these two 
quantities are needed to close Eq. (30). From Eq. (31) one gets 



Equation (33) can alternatively be derived directly from Eq. ( 1 3) by applying the approximations 



d\ I d(0\mf -^ nuii, '- I d\ I dmyco'-f -^ nuT'°\ 



(33) 



(34) 



which are weaker than Eq. (32). 

The simple BGK-like model (30) is formally analogous to the one proposed in Ref. [5] for smooth spheres, Eq. (12). 
The key difference is that the cooling rate ^ is replaced by the energy production rate ^"^, which depends not only on 
r'' but also on £1 and T™', the latter two quantities obeying Eq. (33). Despite its crudeness, this model can be useful 
to explore the basic influence of roughness on the translational properties of a dilute granular gas. 
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FIGURE 1. Surface plots of —P„,/nT^'' (left) and T^^ /T^ (right) as functions of a and /3 in the simple shear flow of a dilute gas 
of inelastic rough hard spheres with a uniform mass distribution (k= |). 





FIGURE 2. 

inelastic roui 



Parametric plots of (a/Av) and {P„/nT^'^) as functions of T^^ /T^^ in the simple shear flow of a dilute gas of 
;h hard spheres. The circles indicate four representative cases. 





FIGURE 3. Surface plots ofR{y,,Vy) = ^""(Vv, Vy)/g{)'(V,-, Vy) for a = 0.8 and jS = - 1 (left) and for a = 0.8, /3 = 0.2, and (f = § 
(right) in the simple shear flow of a dilute gas of inelastic rough hard spheres. 



APPLICATION TO THE SIMPLE SHEAR FLOW 

The simple (or uniform) shear flow is an incompressible flow defined by a linear velocity field u = ayx, a uniform 
density, and a uniform temperature. This paradigmatic state is macroscopically quite simple since only a hydrodynamic 
gradient exists (a = dux/dy) and moreover it is a constant. On the other hand, it is important to remark that the steady 
state resulting from the balance between inelastic cooling and viscous heating is inherently non-Newtonian [13]. In 
the steady state Eq. (30) becomes [6, 13] 

-aVy^r{\) = -A(a,/3)v [^(V) -/^(V)] + ^^ ■ [vr(V)] . (35) 

Taking into account that (5( + u • V = in the steady simple shear flow, Eq. (33) yields il = 0, i^™' = 0. The latter 
condition determines the ratio between the rotational and translational temperatures from Eq. (25). Insertion of this 
ratio into Eq. (24) gives <^"^. The results are 

Interestingly, the temperature ratio is independent of a. It ranges from 7"™'/^"^ = in the perfectly smooth case 
(j3 = — 1) to 7"™'/ rf = 1 in the perfectly rough case (j3 = 1). This temperature ratio is the same as in the homogeneous 
steady state driven by a white-noise thermostat. Multiplying both sides of Eq. (35) by ViVj and integrating over velocity 
one gets a set of algebraic coupled linear equations for the elements Pij — m/iiVy,y,/''(V) of the pressure tensor. 
The structure is the same as in the smooth case [6, 11] and so only the final results are quoted here: 






where <^"^ = <^"^/Av is the scaled energy production rate and we have introduced the anisotropic temperatures Tj"^ = 
Pii/n, where Pa is the ith normal stress. The absolute maximum value of <^"' is 254 (109— 16\/rT) ~ 0.95 and 

corresponds to a = 0, K = | and j3 = 13 - 4\/TT ~ -0.27. Elimination of ^^ between T^/T^ and either a^/P' 
or Pxy/tiT^ yields two nonequilibrium "equations of state" independent of a and j3, namely 

"'^(^) =2(7^] (^"^j' (#) =2r^V'r^)- ^^^^ 

Figure 1 shows the dependence of the ratios —Pxy/nT^'^ and T^J /T^'' on both coefficients of restitution a and j3 for the 
case K" = § (uniform spherical mass distribution). The anisotropic effects induced by the shearing are measured by the 
departure of P^y/tiT^'^ and r,"^/r"^ from and 1, respectively. We can observe that, for a fixed value of j3, those effects 
increase monotonically with decreasing a. On the other hand, for a fixed value of a the influence of roughness is not 
monotonic, the higher deviations from anisotropy taking place for central values of j3 . The dependence on j3 becomes 
less important as the inelasticity increases. In the elastic case (a — 1) and j5 = ±1, Pxy/riT^ — and Ty'/T^' = 1, as 
expected. The two universal equations of state are plotted in Fig. 2. Every triad (a, j3 , fc) is represented by one single 
point in each graph and all the representative points lie on the curves given by Eq. (38). The end points of the curves 
correspond to (a = 0,j3 ~ — 0.27,k'= |) and (a — l,j3 = ±1,)C — arbitrary). 

A practical advantage of kinetic models is the possibility of obtaining explicitly the velocity distribution function. 
Exploiting the analogy with the smooth case [6], the solution to Eq. (35) turns out to be 



f{y) = fdse-^'-'^^"yf^ (ef^"^ (V + asVyx)) . (39) 



This expression is formally analogous to that of an ordinary fluid with a Gaussian thermostat [3]. The marginal 
distribution g^'^{Vx,Vy) = J^a!,dV^f^'^{\) is also given by Eq. (39) except for the replacements |<^"' -^ t,^, f^iy) -^ 
goi^x^Vy) = J^^dVJ^'iY). Figure 3 show the ratio R{Vx,Vy) = ^"'(Vt, Vy)/go (Vx, Vy) for inelastic smooth spheres with 
a = 0.8 and for inelastic rough spheres with a — 0.8, j3 — 0.2, and K" = i. As can be observed, the distortion from the 
local equilibrium distribution is higher in the latter case than in the former. 



CONCLUDING REMARKS 

In this paper a two-level BGK-like description for a dilute granular gas of inelastic rough hard spheres has been 
proposed. At a more fundamental level, the one-body joint distribution function /(v, 0)) is assumed to obey the kinetic 
equation (13) with the replacement (28). The model preserves the Boltzmann collisional integrals (17)-(19) with the 
de-spinning rate (^q) and the partial energy production rates (S,^ and <^™') given by Eqs. (23)-(25). At a simpler and 
cruder level, the translational distribution function /"^(v) is assumed to obey Eq. (30), complemented by Eq. (33), while 
the rotational distribution /"'(o) satisfies Eq. (31). The second level prevents one from accounting for correlations 
between the translational and rotational degrees of freedom [9]. On the other hand, the model kinetic equation (30) 
can be useful to investigate the basic influence of roughness on the translational properties of the granular gas. 

The kinetic model made of Eqs. (30) and (33) has been applied to the steady simple shear flow problem. The solution 
predicts that the temperature ratio 7"™y7""' is independent of the coefficient of normal restitution a and is given by 
Eq. (36). On the other hand, according to Eq. (37), the reduced shear stress Pxy/nT^, the reduced shear rate a/Xv, and 
the anisotropic temperature ratios T^'^/T^'^ and Ty /T^'' depend on a and j3 only through the reduced energy production 

rate <^''. It is observed that, at a given value of a, the dependence of those quantities on j3 is not monotonic, this effect 
being less pronounced as inelasticity increases. 

It is planned to carry out computer simulations to test the above theoretical predictions for the simple shear flow. 
Moreover, the kinetic model (30) will be applied to other states that have been solved in the case of perfectly smooth 
spheres, such as the Couette flow [12], the gravity-driven Poiseuille flow [14], and the uniform longitudinal flow [15]. 

ACKNOWLEDGMENTS 

The author is grateful to V. Garzo and G. M. Kremer for insightful discussions. This work has been supported by the 
Ministerio de Educacion y Ciencia (Spain) through Grant No. FIS2007-60977 (partially financed by FEDER funds) 
and by the Junta de Extremadura (Spain) through Grant No. GR10158. 

REFERENCES 

1. C. Cercignani, The Boltzmann Equation and Its Applications, Springer Verlag, New York, 1988. 

2. P. L. Bhatnagar, E. P. Gross, and M. Krook, Phy. Rev. 94, 511-525 (1954); P Welander, Akiv for Fysik 7, 507-553 (1954). 

3. V. Garzo and A. Santos, Kinetic Theory of Gases in Shear Flows. Nonlinear Transport, Kluwer Academic Publishers, 
Dordrecht, 2003. 

4. N. Brilliantov and T. Poschel, Kinetic Theory of Granular Gases, Oxford U. P., Oxford, 2004; C. S. Campbell, Anna. Rev. 
Fluid Meek. 22, 57-92 (1990); J. W. Dufty, J. Phys.: Condens. Matt. 12 A47-A56 (2000); I. Goldhirsch, Anna. Rev. Fluid 
Meek. 35, 267-293 (2003). 

5. J. J. Brey, J. W. Dufty, and A. Santos, / Stat. Phys. 97, 281-322 (1999). 

6. A. Santos and A. Astillero, Phys. Rev. E 72, 031308 (2005). 

7. S. J. Moon, J. B. Swift, and H. L. Swinney, Phys. Rev E 69, 031301 (2004). 

8. A. Zippelius, Physica A 369, 143-158 (2006). 

9. W. T. Kranz, N. V. Brilliantov, T. Poschel, and A. Zippelius, Eur Phys. J.- Spec. Top. 179, 91-111 (2009). 

10. A. Santos, G. M. Kremer, and V. Garzo, Prog. Theor Phys. Suppl. 184, 31-48 (2010). 

11. J. J. Brey, M. J. Ruiz-Montero, and F. Moreno, Phys. Rev E 55, 2846-2856 (1997). 

12. M. Tij, E. E. Tahiri, J. M. Montanero, V. Garzo, A. Santos, and J. W. Dufty, / Stat. Phys. 103, 1035-1068 (2001). 

13. A. Santos, V. Garzo, and J. W. Dufty, Phys. Rev E 69, 061303 (2004). 

14. M. Tij and A. Santos, / Stat. Phys. 117, 901-928 (2004). 

15. A. Santos, "Longitudinal Viscous Flow in Granular Gases," in Rarefied Gas Dynamics: Proceedings of the 26th International 
Symposium on Rarefied Gas Dynamics, edited by T. Abe, AIP Conference Proceedings, vol. 1084, Melville, NY, 2009, pp. 
93-98. 



